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Abstract 

A "bubble universe" nucleating in an eternally inflating false vacuum will experience, in the 
course of its expansion, collisions with an infinite number of other bubbles. In an idealized model, 
we calculate the rate of collisions around an observer inside a given reference bubble. We show 
that the collision rate violates both the homogeneity and the isotropy of the bubble universe. Each 
bubble has a center which can be related to "the beginning of inflation" in the parent false vacuum, 
and any observer not at the center will see an anisotropic bubble collision rate that peaks in the 
outward direction. Surprisingly, this memory of the onset of inflation persists no matter how much 
time elapses before the nucleation of the reference bubble. 
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I. INTRODUCTION 

In a theory with different metastable vacua, such as the landscape of string theory 
the process of eternal inflation may lead to a "multiverse" where different regions of space- 
time are occupied by different vacua. A region occupied by vacuum A may spawn regions of 
other vacua B,C,D... which are adjacent in field space. Daughter regions will in turn spawn 
their offspring, and so on, producing an infinite tree of "pocket universes." 

Vacua can evolve by a mix of semiclassical tunneling and stochastic evolution, but in this 
paper we will consider only tunneling. In that case, bubbles form by random nucleation and 
then start to expand into the parent vacuum with constant acceleration. The interior of the 
growing bubble has the geometry of an open FRW universe |3] , and (assuming that a short 
period of slow roll inflation flattens it out to satisfaction) we may entertain the possibility 
that we live in one of such "bubble universes" {4]. 

The above description, however, is incomplete because it ignores collisions with other 
bubbles. Collisions may be quite rare if the nucleation rates are small. Nonetheless, a 
bubble expands for an infinite amount of time, and will collide and merge with an infinite 
number of other bubbles, forming an ever growing "cluster" . 

We should therefore reevaluate the naive picture of a smooth FRW universe on large 
scales. In particular, we may ask what fraction of the FRW universe remains unaffected by 
collisions. As we shall see, only a set of measure zero remains unaffected. If this is supposed 
to describe our local universe, there seems to be reason for concern. Why haven't we seen 
any collisions yet? Do we have much time left until we are blown away by a collision? What 
is the expected distance to the nearest point in our FRW time slice which has already been 
hit by another bubble? The purpose of the present paper is to explore some of these issues. 

In the course of this investigation we have stumbled upon a rather remarkable result. We 
find that the rate of collisions around a typical observer in the bubble universe is anisotropic, 
and the origin of this anisotropy is related to the beginning of false vacuum inflation. 

A metastable inflationary de Sitter phase can be eternal to the future, but not to the 
past. More precisely, the inflating region of spacetime is geodesically past-incomplete ij, 
and therefore some initial conditions must be specified on the past boundary of this region. 
For instance, we may posit that at some initial time, a given large region of space is in false 
vacuum. The congruence of geodesies normal to the initial time surface defines a preferred 
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frame, with respect to which velocities can be defined. 

The standard lore is that this preferred frame is not important, and that the initial 
conditions are soon forgotten. This may be true for geodesic observers in the false vacuum 
phase, whose velocity with respect to the preferred congruence redshifts exponentially with 
time. For those "co-moving" observers, memory of the initial surface is lost as we push the 
surface far away into the remote past. However, when a bubble of a new phase forms in 
the original false vacuum, the congruence of observers in the FRW open universe includes 
observers with velocities arbitrarily close to the speed of light relative to the preferred 
congruence. Such individuals are not "far" from some points on the initial surface, and may 
have a chance to detect some information about the beginning of inflation. 

The plan of the paper is the following. In Section II we describe the geometry and statis- 
tics of collisions onto a given reference bubble. In Section III, we consider the anisotropics 
in the bubble distribution as seen by observers inside a reference bubble. The origin of these 
anisotropies is further discussed in Section IV in a simplified context, where the observer is 
in false vacuum rather than inside the bubble. Section V is devoted to conclusions. Some 
technical details are left for the appendices. 



II. COLLISIONS ONTO A REFERENCE BUBBLE 

Here, we shall derive some basic results concerning the distribution of collisions impinging 
on a given reference bubble. We begin with a description of the basic setup we shall consider. 

It is convenient to use flat de Sitter coordinates to describe the background inflating false 
vacuum, 

ds 2 = de-e 2Ht {dr l + r l dto 2 ). (1) 
To simplify the equations, we shall choose units so that 

H = l. (2) 

To simplify matters further, we shall assume that the vacuum energy density inside the 
bubbles is nearly the same as that outside (at least for a sufficiently long time after bubble 
nucleation) and that the gravitational effect of bubble walls is negligible. Then the metric 
is 

ds 2 = dt 2 -e 2t (dr 2 + r 2 dQ 2 ) (3) 
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in the entire spacetime region of interest. 

It is also useful to consider the "embedding" of de Sitter space as a timelike hyperboloid of 
unit radius, in a 5 dimensional Minkowski space, whose rectangular coordinates are labeled 
V,W and X = (X,Y,Z): 

X 2 + W 2 -V 2 = 1. (4) 
These are related to the flat chart coordinates t and x = (x,y,z), (r = \x\), through 

W - V = e - * - eV 2 , W + V = e\ X = e t x. (5) 

Let us consider a "reference" bubble that nucleates at t = r = (we shall sometimes call 
this "our" bubble). In the embedding coordinates, this corresponds to 

W = l, V = X = 0. (6) 

The bubble geometry is symmetric under boosts which have the nucleation event as a fixed 
point. These form an 0(3, 1) group of isometries. 

The interior of the bubble (or more precisely, the interior of the light cone from the 
nucleation event) is described by the line element 

ds 2 = dr 2 - sinh 2 r (d£ 2 + sinh 2 £dtt 2 ) . (7) 

The coordinates (t, r) and (r, £) are related by 

e* = cosh r + cosh £ sinh r, (8) 
eV = sinh £ sinh r. (9) 

We set the initial condition that there are no bubbles at some t — U. This breaks the residual 
0(3, 1) invariance of the bubble, which is responsible for the homogeneity and isotropy of 
our FRW universe. Consequently, not all observers who live in the open FRW universe will 
see the same. 

A. Collisions around the observer at £ = 0. 

Let us now concentrate on the distribution of collisions around the point £ = £ Q b s = 0, 
which is at the origin of the open FRW hyperboloid. We leave the consideration of typical 
observers, who live far away from the origin, to the next Subsection and to Section III. 
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The distribution around £ Q b s = will of course be isotropic. In particular, we shall be 
interested in the typical distance at which we might expect the nearest collision. In what 
follows, we shall calculate the probability P(£, r) that no collisions with other bubbles have 
affected a spherical region of radius £ around the origin, on a hypersurface of constant r. 

As a warm up exercise, let us consider the following question. Assuming that the point 

t = r = is still in false vacuum, what is the probability that some bubble will have hit 

the surface t — at some r < ro?. The relevant quantity is the 4-volume inside the past 

light-cone of the the circle r < r minus the 4-volume inside the past light cone of the 

origin r = (see Fig. 1). The physical radius of the past light cone from the origin is 

Ro(t) = 1 — e* (t < 0), whereas the radius of the past light cone from the circle of radius r 

is Ri(t) = Ro(t) + r e*. The volume of the grey shaded region in Fig. [T]is thus given by 

47T f° 2tt 
V 4 (r 0> ti) = y / OR? - R 3 ) dt = — (6r + 3r 2 + 2r 3 ) + O(eV ). (10) 

Note that this is finite (for finite ro) even in the limit when the initial surface is pushed all 

the way to — > — oo. For a region of a Hubble size, we have 

. 22vr 
V 4 (l,tj -> -oo) 
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Hence, the probability of having a bubble one Hubble distance away from ours at t = is 
of order of the nucleation rate per unit volume A, which we shall assume to be small (A < 1 
in the units where H = 1). The distance to the nearest bubble at t — can be estimated 
from the condition AV4 ~ 1. From Eq. (TlPl) . this distance will be of the order ro ~ A -1 / 3 . 

Next, we may consider the same question on a r = const, hyperboloid. Past directed 
radial null rays define a mapping between points on the hyperboloid at some distance £ and 
points on the "bubble cone," by which we mean the future light-cone from the nucleation 
event. The bubble cone is given by the equation 

r(t) = l-e -t . (t>0). (11) 

The flat chart coordinates of points on the hyperboloid can be found from Eq. (Q. The past 
directed outward radial null geodesic from a point (£, r) is given by 

r(t) = e- t + F(r,Z), (12) 

where 

sinh £ sinh r — 1 , . 

(13) 



cosh r + cosh £ sinh r 
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FIG. 1: Region in grey shade, V^q) can nucleate bubbles which will be within the distance ro 
from the origin (on the flat surface t = 0). Here, we assume that the point r = t = is still in 
false vacuum. Region in horizontal stripes, V4(£,t), can nucleate bubbles which will collide with 
ours, and which will be visible within a coordinate radius £ from the origin (on the r = const. 
hypersurface) . Here, we assume that the center point £ = has not yet been hit by any bubble. 

Note in particular that for £ = 0, we have F = —e~ T , whereas for large £ and r, we have 

F{£, t) w 1 - 2e"* + 2e~ 2€ - Ae'^. (£,r»l) (14) 



The intersection of the null geodesic (1121) with the bubble cone ( TTT1) is at 

e - = iz£, r ..l±£. (15) 
(Incidentally, the intersection with the plane at t = is at r = 1 + F < 2.) 
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It will be very convenient to use Eq. ffT2l) as the definition for a change of variables, 
to replace r in favor of F. The different values of F can be thought of as labeling the 
different past directed outward radial null rays emanating from a given £, on a r = const. 
hyperboloid, as indicated by Eq. (fT5|) . More explicitly, we introduce 

a = e\ F = r-e~\ (16) 

in terms of which the metric reads 

ds 2 = IdadF - a 2 dF 2 - (1 + aFfdtt 2 . (17) 

This form of the metric is completely regular at a = (which corresponds to the boundary 
of the flat coordinatization). In fact, the new chart (j!7p covers the whole of de Sitter space, 
when we let the coordinates vary in the range — oo<a<oo, — oo<F<oo, with the 
restriction aF > —1. Negative values of a correspond to the part of de Sitter space not 
covered in a flat chart. 

Let us assume that, on a given r = const, hyperboloid, the origin £ = has not yet 
been hit by a bubble, and let us ask what is the typical distance at which we may expect 
the nearest collision. The probability that no collisions have affected a region of coordinate 
radius £ around the origin is given by 

P(£,T,t i ) = e- xv \ (18) 

where A is the nucleation rate, and V4 is the relevant 4 volume to the past of a spherical 
region of radius £ on the hyperboloid. Here, we must subtract the contribution from the 
interior of the past light cone of £ = (this is because we assume £ = has not been hit 
by any bubbles), and also the intersection with the interior of our bubble, since we assume 
that no new bubbles can nucleate inside of our bubble. We may express this 4 volume as 

*( { , r, t( ) = 4, / ^« = 4, / e=V< + F) 2 dFdt = 4, /(I + aFf^F. (19) 

The evaluation of V4 is simplest in the (F, a) coordinates. This is because the light cones 
are bound by null geodesies with F = const. Hence, the integration limits for F in ( Tl9l) are 
independent of a, 

-e~ T < F < F(r,0- 
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The lower limit corresponds to the center of the hyperboloid £ = 0. The integration range 
for t depends on F: £j <t < t c (F), where t c is given in ffl5l) . Here we are approximating the 
boundary of our bubble by the "bubble cone," which is justified when the size of bubbles at 
the time of nucleation is much less than the Hubble radius. In terms of a — e*, the range is 
given by 

Oi < a < 2/(1 - F). 
Performing the double integral, we have 



2F 



ln(l — F) + 0{ ai F) 



(20) 

F=-e~ T 



l-F) 

where the upper limit is given by (fT3"l) . Note that the 4- volume takes a finite value in the 
limit when the initial surface is pushed into the remote past tij — ► 0. 
Defining K 4 (£, r) = lim^.oo 1/ 4 (£, r, t { ), we find 



Vi(£, r) = y je 5 sinh£ tanh 2 + In 



1+^(1-0(^-1; 



(21) 



In particular, for r ^> 1, we have 



V^t) « 4< + 0(a (£«1), (22) 
V^,t) « y e 2 « + O(0 (£»1). (23) 

The typical coordinate distance to the nearest bubble is found from AV4 ~ 1, and for small 
nucleation rates, A -C 1, it is of order 

A£~ln(l/A). (24) 

Note that this is independent of r. This is to be expected, since (for r > 1) the products 
of collision are "frozen in" with the expansion of the universe during the period of inflation 
inside the bubble. Here, for simplicity, we have approximated this period as a de Sitter 
phase with the same expansion rate as the false vacuum. In a more realistic case we should 
consider a period of slow roll inflation, followed by a standard decelerating phase. If this 
phase leads to a Minkowski vacuum, all products of collision would gradually fall within the 
horizon, and eventually hit any observer inside the bubble. (This is explicitly shown at the 
end of Appendix 1, where we calculate the rate at which an observer in Minkowski vacuum 
would see new bubbles falling into his field of view.) On the other hand, our own universe 
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does not seem to be approaching a Minkowski vacuum. Rather, it appears to be accelerating 
again, and collisions with other bubbles may forever remain hidden behind the cosmological 
horizon. 

The curvature radius of the FRW universe corresponds to A£ ~ 1. Current observational 
bounds on the spatial curvature imply [6] (a H )~ 2 = |1 — fi| < 10~ 2 , where the subindex 
indicates the present time. The coordinate size of the observable universe is given by 
A£ ~ (cto-f^o) -1 - Hence, the observer at the center of the hyperboloid £ bs = can only 
see out to coordinate distances A£ < 10 _1 . Using ([18]) and (|22|) . the probability that this 
observer may see any collisions at all is given approximately by 

1 - P(A£ ) « 1 - e- 47rAA?0 < A < 1. 

This issue will be discussed further in Section III, and also in Appendix 1. In the Appendix 
we show that the asymptotic expressions (1221) and (|23|) are valid also for the relevant four 
volume in the neighborhood of any observer (at £ Q b s 3> 1) who has not yet been hit by a 
bubble. This means that we are very unlikely to be hit by a bubble at any time in the 
future. 

The expected physical distance to the collision nearest to £ = 0, at the time when inflation 
ends, can be estimated as 

d = Zln(l/A), (25) 

where Z ~ sinhr/ is the slow-roll expansion factor (ry is the time at the end of inflation 
inside the bubble). Here we have worked in a simplified scenario where H is the same in 
false vacuum as it is in the period of slow-roll inflation inside the bubble. In the general 
case, the expression for d will be more complicated, but the dependence on Z and A will be 
similar. 

B. Fractal dimension of the bubble universe. 

The world-line of a given point (8, <fi) on the wall of the reference bubble will sooner or 
later be hit by other bubbles, with probability equal to one. This follows from the fact that 
there is a finite probability per unit time t for this worldline to be hit by other bubbles (see 
Section HVl for a rigorous derivation of this statement). Hence, the fraction of the bubble 
wall area that makes it to future infinity without collisions is a set of measure zero. 
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The effect of collisions onto the reference bubble propagates into the open FRW universe. 
On the hyperbolic slice r = const., each collision affects a "wedge" shaped region, whose 
tip points towards the origin and which extends all the way to £ — > oo, spanning a finite 
asymptotic solid angle at spatial infinity Let us denote by £ c the distance from the tip of 
the wedge to the origin, and let us calculate the asymptotic angular size of the wedge, 9 W , 
as a function of £ c . 

It is convenient, again, to use the coordinates a, F of the chart (fTTj) . A bubble that nucle- 
ates at time a has asymptotic co-moving radius 1/a at future infinity. A bubble nucleating 
at point F will have its center displaced a distance r = F + a~ l from the origin, and its 
boundary at future infinity will have co-moving cartesian coordinates satisfying: 

x 2 + y 2 + (z - F - a- 1 ) 2 = a~ 2 . (26) 

(Here we have assumed that the center is on the z axis). The future boundary of the 
reference bubble (which nucleates at t = r = 0) is given by 

x 2 + y 2 + z 2 = 1. (27) 

The intersection of both is at 

la + 2F + aF 2 

z = cos9 w = (28) 

2 1 + aF K ' 

Note that 9 W corresponds to the asymptotic angular size of wedges corresponding to bubbles 

nucleated at (a,F). 

On a hyperboloid with r ^ 3, large £ corresponds to 1 — F w 2e~^(l + coshr)/ sinhr pa 
2e~^. From (TT9l) . the bubbles which affect this region nucleate in a four- volume dominated 
by the largest possible a, 
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From (Izol) . The typical angular size of wedges whose tip is at £ c 3> 1 is therefore 

e w (Q ~ « 1. (29) 

Let us denote by i\(£) the solid angle which remains unaffected by collisions out to a 
distance £ from the origin. As £ is increased, more and more wedges accumulate, each one 
removing a solid angle 2ir[l — cos^] from Q u . The loss of solid angle due to wedges whose 
tip is in the interval d£ is given by 

dSl u (£) ^ -2tt[1 - cos9 w ]XdV 4 ~ -Xe^dV A , (30) 
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where dV^ is the region of four volume to the past of the interval g?£, in the unaffected 
portion of the sphere. Here, we are making the approximation that wedges do not overlap 
with each other. Also, we shall assume that a wedge whose tip is at £ c depletes a solid angle 
equal to the corresponding asymptotic value (at £ — > 00) for all £ > £ c . The four volume is 
approximately given by Eq. (j!9p . where 4ir is replaced by i\(£). Using (j28j) . we have 

-2/(l~F) 



2tt[1- cos 6 w }dV 4 = 7i-n u (0 



/ (1 - F)[2 - a(l - F)](l + aF)da 

Ja=0 



dF 



= ^_^ n u(0 dF « fn u m, (31) 

where in the last step we have used F ~ 1 — 2e~£. 

Combining Eq. (I3TT) with (1301) . the solid angle which remains unaffected by collisions as 
a function of £ is given by 

fi u (0~4vre-^, (£»1), (32) 

where k ~ (47r/3)A. The volume element on a spatial section of an open FRW is given by 
dV FRW = 47rsinh 2 £d£. The volume which is unaffected by bubbles 

= fi u (£) sinh 2 £d£ ~ 27re (2 - K) «rf£, (33) 

is therefore unbounded and dominated by large distances from the origin (assuming, of 
course, that /«<1). Nevertheless, the unaffected volume fraction tends to zero in the limit 
£ — > 00. The unaffected volume fraction will be recalculated in the next section, using a 
method which is exact in the context of our idealized model. 

The fact that the unaffected part of the volume is a set of measure zero might seem to be 
reason for concern. Nevertheless, if we pick a random point in the unaffected region, then 
the closest hit is likely to be quite far away. Indeed, the number of bubbles hitting in the 
interval g?£ is given by dN = XdV^. Using (T211 . fl3T|) and fl3"31 we find that for large £, 

dN ~ XdV. (34) 

Let us now pick a random point in the unaffected region, at £ 3> 1, and let us choose it as 
the origin of open FRW coordinates. A spherical region of radius A£ ^> 1 around that point 
has coordinate volume V ~ e 2A ^. Using (13~41) . the expected distance to the nearest bubble 
will be of order 

A£ ~ ln(l/A). (35) 
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This is the same result we found in the previous Subsection, for the distribution around 
the privileged point at the center of the hyperboloid. The coordinate distance fl35|) is large 
enough that we should not be too concerned about the hazard of future collisions. 

Finally, let us characterize the fractal dimension of the unaffected part of the bubble uni- 
verse. Looking outward from the center of the hyperboloid, we see finer and finer "wedges" 
carving away the solid angle as we increase £. If we define a smearing angle e = e - ^, and 
ignore structures smaller than e, the unaffected solid angle is fi u (£) ~ 47re K . The number of 
sets of angular radius e which is needed to cover this region is n ~ e K ~ 2 . Thus we can think 
of 
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D = 2- k = 2 —A (36) 

as the fractal dimension of the unaffected portion of the "celestial" sphere at £ — > oo. Clearly, 
this is also the fractal dimension of the remaining surface of our bubble which has not been 
hit by other bubbles at future infinity. 



III. ANISOTROPIES IN THE DISTRIBUTION OF COLLISIONS 

To study collisions around points which are away from the origin, it is useful to perform 
a de Sitter transformation (Lorentz transformation in the embedding space) to a new frame 
S' where the point of interest is at the origin. This greatly simplifies the geometry of the 
relevant past light cones. 

If the original point at large £ = £ D b s was along the Z direction, then we use the boost 

V' = j(V-PZ), Z' =j(Z-PV), X' = X, Y' = Y, W' = W } (37) 

with (3 = tanh£ b s and 7 = cosh£ b s , in order to bring the point of interest down to £ obs = 0. 
In doing so, the initial surface t = ti gets distorted and this is of course something we have 
to consider. 

In the embedding coordinates, t = ti corresponds to the null plane 

W + V = e u = Oi, (38) 

which in boosted coordinates reads 

V' + PZ' = 1 -\a i -W). (39) 
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We can now express this in terms of flat chart coordinates (t', x') by using the standard 
relations (jSJ). This leads to 

sinht' + 2 e< ' r ' 2 + /3eV cos0' = 7 _1 (a* - cosht' + ^eV 2 ). (40) 

In terms of coordinates (a', F') analogous to the (a, F) pair which we introduced in the 
previous Section [see Eq. (JT71) ]. the equation for the initial surface a = a\ takes the more 
tractable form a! = a'^F', 9'; cti), where 

j 2 q,-/3 7 cosg'-( 7 -l)F' 

l + 7 + 2/3 7 F'cos^ + (7-l)F' 2 ' 1 ' 

It should be noted that a[ can be negative for some values of F' . This is not a problem. 
It just means that in the boosted frame, the initial surface invades the portion of de Sitter 
space not covered in the original chart. 

It is interesting to calculate the probability that the FRW observer at £' = (that is, 
£ = £obs) will be affected by a collision with another bubble before some specified time r. 
This is related to the spacetime volume available for the nucleation of bubbles which lies in 
the past light cone of the observation point (£' = 0,r), but which does not lie inside our 
reference bubble, or inside the past light cone of the nucleation event. This four volume 
increases with proper time, in the following way 

" dF' 



dV* 



2/(1-*") 

(1 + a'F'fda! 

a' 



dr 



■drdQ', (42) 



where dQ! = 2n d(cos9'). The probability per unit proper time and solid angle that our 
observer is hit by a bubble, assuming that he was not already hit by a bubble, is given by 

dP , dV± 

^^kr (43) 



dr dVL' drdVt' 

From Eq. (1131) . £' = corresponds to F' = —e~ T , and performing the integration over a' we 
find 



dV* 1 



drdQ' 3 



cosh £ obs sinh r + cosh r — Oj ^ 1 3 ( r 



- tanh d 



(44) 



cosh £ b s cosh t + sinh r — sinh £ Q b s cos 9' J V 2 , 

Note that the result depends both on the point of observation £ b s and on the direction of 
observation 9'. Homogeneity and isotropy are lost, even in the limit when we push the initial 
surface all the way to U — > — oo, i.e. eij = 0. 
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For any fixed nonzero value of r, the limit of (|Hj) as £ Q b s — » oo is given by 



sinh r 



cosh r — cos 9' 



tanh 3 ( - 



(45) 



The existence of this limit is good news. The number of observers grows without bound with 
the distance to the center of the hyperboloid, and thus we expect that the typical observer 
lives at very large £ b s - He or she should therefore measure the distribution (145]) . Also, the 
limit is independent of £j. This is also interesting, since it means that some specific details 
about the surface of initial conditions do not seem to matter. Nevertheless, the distribution 
is anisotropic, with the minimum number of hits per unit time in the direction of the center 
of our reference bubble. In this sense, memory of the initial surface persists. For large r, 
the distribution takes the simple dipole form 

dV A 



drdtt' 



2(l + cosfl')e~ r + 0{e- 2r ). 



(46) 



Dipole anisotropics and memory of initial conditions are not commodities one usually 
expects from inflation. Although we find this result to be rather shocking, the effect is real. 
Its origin is best understood by eliminating the complications due to the bubble geometry, 
as will be discussed in the next section. First, however, we would like to study a bit further 
the implications of ( 1441 ). 

To find the total rate XdV^/dr at which bubbles will be encountered by an observer who 
has not previously been hit by a bubble, one can integrate (jHJ) over solid angle. The result 
is given by 



dV 4 

~d7 



Air 

T 



(cosh £ bs sinh r + cosh r — a^) 3 (cosh £ Q b s cosh r + sinh r) 



- tanh 3 



T 



(47) 



(cosh £ obs sinh r + cosh r) 4 V 2 , 

The dependence of this result on £ obs shows us the inhomogeneity of the bubble collision 
rate. For the special case of r = the expression simplifies to 

dV A An 



-3— = — (1 - (k) cosh£ obs , 
dr 6 

which is highly inhomogeneous, while for large r the expression becomes 

dV A „ ai 



57T 



dr \ cosh £ obs + 1 

which shows that the inhomogeneity disappears when t 
when £ obs is large. 



e- T + 0(e- 2T ) 



-oo i.e., a,- 



(48) 
(49) 

0), and also 
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We can continue by integrating (j4T)l over r, from zero up to an arbitrary value, thereby 
determining the total 4-volume available for the nucleation of bubbles that could collide 
with an observer located at (£ b s ,T). One finds 



V 4 



4tt 



< tanh 2 



3 1 V2. 

- 1~ < cud - r ) 



In 



e cosh 



where 



-aj 1 



1 



(50) 



(51) 



cosh £ b s sinh r + cosh r 
This quantity represents the 4-volume of the region that is in the past light-cone of the point 
(£obs) T )) but is not inside the reference bubble nor in the past light-cone of its nucleation 
event. This volume can be calculated directly without making the coordinate transformation 
(I37p . and we have verified that the results agree. 

If an observer at £ = Us has not seen a bubble at time T\, then the probability that she 
will not be hit by a bubble by some later time r 2 is determined by (loDj) . with 



P = exp {-A [V^fobs, t 2 ) - V^obs, n)}} . 



(52) 



It is particularly interesting to look at IjoOj) for large Us, since a typical observer will be 
found at arbitrarily large Us- In that limit one finds that 

2 e -?obs 4 cQgh T 



C — — o 

smh r sinh r 
which is valid for any r^O. Then 

V 4 = y jUs + tanh 2 (^j + In 



O (e 



tanh ( — 
2 



1 2 1 o e" 
-47T < a* a,- H — aj 1 — 2aj- 



2 1 ' 9 % 'sinhr 
Then, if we are interested in the limit r —>■ oo, we find 



+ 2e 5obs coth(r) 
+ O(e- 2?obs ). 



lim ^ 4 = % (Us + 1 + 2e~^) - 4tt f Oi - \a\ + i a H + 0{e 
r^oo 3 \ 2 9 / 



2?obs 



(53) 



(54) 



(55) 



The probability of a point not being hit by a bubble is given by P = exp(— AV4), so the 
leading term 47r£ bs/3 in the formula above reproduces Eq. ( 1321) . which was used to determine 
the fractal dimension. 
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IV. A MOVING OBSERVER IN DE SITTER 



Even for an observer in the false vacuum, the bubble nucleation rate and the angular 
distribution of bubbles depend on the observer's velocity relative to the "preferred" co- 
moving congruence C which is determined by the surface of initial conditions. Here, we give 
a self-contained account of this effect. 

Consider a de Sitter space with H — 1, 

ds 2 = r]- 2 (dr] 2 -rfx 2 ), (56) 

with — oo < i] < 0. The conformal time rj is related to the usual time vatiable as 

r] = -e~\ (57) 

We shall assume that there are no bubbles at some initial moment rj = rji. 
Consider an observer at x = 0, r\ — t]q, moving with a velocity 

v = tanh0 (58) 

relative to the comoving observers of (I56I) . We want to know the probability for this observer 
to be hit by a bubble per unit proper time (by his clock). 

An infinitesimal proper time interval St corresponds to a coordinate displacement 

Si] = \rjo\ST cosh 0, (59) 

5x = vSrj. (60) 

The probability to be hit by a bubble is determined by the spacetime volume between the 
past light cones of the points (rj , 0) and (r] + Sr], 5x). The first of these light cones is given 
by 

|xi(77)| = Vo -V (61) 

and the second is given by 

|x 2 - <5x| = Si] + 770 - 77. (62) 

To linear order in St, 

|x — <5x| ~ r — I rjo \ St sinh <ft cos 9 = r — vSrjcosO, (63) 
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where r = |x| and 9 is the angle between x and v, so we can rewrite (|62|) as 

^2(^7, 9) = 5r](l + v cos 9) + r] — 77. (64) 
The spacetime volume between the two light cones is given by the integral 



SV a = 2tt dr]r]- 4 d6 sin 9r{ (r])5r(9) = 2ttSt] / drfrf^fo-ri) l \ d0sin0(l + i;cos0), 

JO JO 

(65) 

where 5r(0) = ^(77, 0) — ri(^) = 5r/(l + vcos9). The integration over 77 is easily done by a 
change of variable £ = — 1/77, 

rvo r 1 

/ cir/r/- 4 ^ - r/) 2 = / + l) 2 = ^0 - 6) 3 - (66) 

Substituting this into (1651) and using Eq. (|59|) for £77, we have 

27T Z* 71 " 

5V 4 = — Sr cosh (pf (t -U) \ d0sin0(l + ucos0), (67) 
3 Jo 

where 

/(*) = (l-e-*) 3 (68) 

and t is the usual time variable which is related to 7/ as in Eq. (|57jl . 

Note that /(t) — > 1 as i — > 00. In what follows we assume the limit ti — > — 00 and set 
f(to - U) = 1. 

The 0-integration in fl67|) is, of course, easily done; the result is 

SV 4 = ^-St, (69) 

where 

St = l^ol" 1 ^ = St cosh (70) 

is the interval of time t corresponding to the proper time interval St. 

The result (1691) is easy to understand. We have two past light cones, one has its origin at 
(to, 0) and the other has its origin shifted in both time and space directions. Eq. f)69p tells 
us that the 4-volume difference SV4 depends only on the time displacement St. And this is 
as it should be. The volumes of the light cones do not change when we shift them in the 
"horizontal" (space) directions, so SV4 remains unchanged, as long as one light cone remains 
entirely within the other. 
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All observers see the same bubble nucleation rate per unit time t, but this corresponds 
to different rates per unit proper time r, 

^ = y cosn 0- (71) 

There is thus a preferred frame, where = 0, corresponding to the lowest nucleation rate. 

The distribution of the arrival directions of the bubbles is also anisotropic. The angular 
distribution can be easily read from Eq.( l6Tj) . 

^^ = l C osh0(l + t;cos#). (72) 

The velocity of a geodesic observer relative to the comoving frame decays as 

v(t) oc e~* (73) 

at large t, so the anisotropy rapidly disappears and the rate approaches that for a comoving 
observer. However, when we study bubble collisions and construct FRW coordinates inside 
a bubble, the set of comoving observers associated with these coordinates includes observers 
with arbitrarily large values of v. This is the origin of the anisotropy measured by typical 
observers in a bubble universe. 

The frame-dependence of the bubble arrival rate is perhaps not too surprising. To the 
moving observer the initial surface t = ti — > — oo looks rather odd. It is a null surface which 
crossed the worldline of the observer a finite proper time r in the past. This time has been 
calculated in [5[; it is given by 

1. (l + l 
r — - In 



2 \7-l. 

For (7 — 1) <C 1, the initial surface is many Hubble times ago, and the difference from a 
comoving observer is small. However, for 7 ^> 1, the cutoff surface is only a small fraction 
of the Hubble time away. No wonder it changes the rate and introduces an asymmetry. 

The connection between (1721 and the results of the previous Section is not completely 
straightforward, and deserves some comment. The main difference is that observers inside 
the bubble will not have any bubbles nucleating in their immediate neighborhood (since we 
assume that bubbles can only nucleate in false vacuum). 

Near the surface r = 0, corresponding to the origin of time in the open chart, this 
difference becomes irrelevant. For r = 0, Eq. (j44p gives 

(74) 



dV 4 1 



dTdW 3 7 3 (1 -pcosO'Y 
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At first sight, this looks rather different from (1721) . However, if we take into account the 
aberration of the angles in a moving frame, 

cost) = — , 

l-/3cos0'' 

we find 

dVi i 7 (l + / 3cos#). 



drdVt 3 

This agrees with (1721) under the identifications 7 = cosh</> and f3 = v. From Eq. (jSJ), the 
FRW observer at constant £ has a relativistic factor with respect to the preferred "rest 
frame" congruence C which is given by 

dt sinh t + cosh £ cosh r 

7 = — = . (75) 

dr cosh t + cosh £ sinh r 

For r = 0, we have 7 = cosh£, which is unbounded for large £. Hence, right after nucleation, 
observers in this FRW congruence have arbitrarily large relativistic factors. 

For any r > 0, the limit of very large £ gives a finite 7 ~ cothr. This means that a 
period of inflation inside the bubble slows down all observers to non-relativistic speed (for 
t ^> 1, the velocity is of order e~ T ). Nevertheless, to the observers at large £, the initial 
surface looks "slanted," sloping down in time in the direction away from the origin. This 
leads to the anisotropy in the distribution of bubbles. In embedding coordinates, the initial 
surface at t — ti, is described as the intersection of the null surface W + V = e* j with the 
de Sitter hyperboloid (jlj) [the coordinate V is "time," whereas the W axis is normal to 
the hyperboloid at the nucleation point (V = 0, W = 1)]. Now, we can always bring an 
observer from very large £ to the origin of coordinates £' = by means of a large boost, with 
(3 = tanh£ w 1. If in the original "rest" frame the observer was at large z, the surface of 
initial conditions in the new reference frame takes the form V' + Z' « 0. This "initial" plane 
is (almost) tangent to the bubble cone, along its null generator in the negative z' direction. 
Since this leaves no room between the initial surface and the reference bubble, it is clear 
that the probability of being hit from that particular direction vanishes. From the point of 
view of the rest frame C the anisotropy is not surprising either. The observers at large £ 
are very close to the light-cone, near the surface of the reference bubble. Other bubbles can 
only hit from the false vacuum outside, and because of that FRW observers are more likely 
to be hit by bubbles approaching from even larger £ than from any other direction. 
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V. CONCLUSIONS 



Our primary goal in this project was to study the effect of bubble collisions on the 
structure of "bubble universes" in the inflating false vacuum. In particular, we wanted to 
know how likely it is for an observer living in one of the bubbles to be affected by such a 
collision. 

In the absence of collisions, the bubble interior is described by an open FRW model. A 
constant-FRW-time slice of such an unperturbed bubble universe is a hyperboloid, a space 
of constant negative curvature. Each bubble collision carves an infinite wedge-like region 
out of the hyperboloid. We found that the part of the hyperboloid that remains unaffected 
by collisions has a fractal character. Its volume is infinite, but it constitutes a vanishing 
fraction of the total volume of the hyperboloid. 

If we pick a random observer in the unaffected region, then we find that the typical 
distance d from this observer to the nearest bubble collision is given by 

d~i21n(l/A). (76) 

Here R is the curvature radius of the hyperboloid and A 1 is the bubble nucleation rate 
(per Hubble volume per Hubble time in the false vacuum). The origin of the logarithmic 
dependence on A is that the volume grows exponentially with distance on a hyperboloid. 
Current observational bounds on the spatial curvature imply R > IOHq 1 , where Hq 1 is the 
present Hubble radius, and Eq. (f76|) yields d > IOHq 1 . 

Assuming that we live in a bubble universe, we have also estimated the probability for a 
collision to occur within our observable range. This is given by 

P cM ~ AuX/HqR < A. (77) 

The bubble nucleation rate A is usually exponentially suppressed, and thus the chance for 
us to observe a bubble collision is rather remote. 

In the process of this investigation, we have uncovered a remarkable fact, that the proba- 
bility for an observer to be hit by a bubble has a strong dependence on the arrival direction 
of the bubble. The origin of this effect can be traced to a simpler setting, which does not 
involve bubble collisions. Consider a geodesic observer who lives in false vacuum. We want 
to know the probability for this observer to be hit by a bubble, per unit time by his clock. 



20 



We found that this probability depends on the observer's velocity v relative to a certain 
preferred "co-moving" frame, which is set by the initial conditions at the beginning of in- 
flation. The bubble nucleation rate is minimal for co- moving observers with v = 0, and the 
most probable arrival direction for a bubble is that opposite to v. 

For a single observer, these effects will be present only for a brief period of time. If the 
initial state of false vacuum is specified on some spacelike hypersurface, the co-moving frame 
is defined by the congruence C of geodesies orthogonal to that surface. The velocity of a 
geodesic observer relative to C redshifts exponentially with time, so the bubble nucleation 
rate rapidly becomes isotropic and approaches its co-moving value. This is in accord with 
the widespread belief that the initial conditions at the onset of inflation have no lasting 
effect. 

We found, however, that this folk wisdom does not apply when the situation is described in 
terms of a FRW open universe inside the bubble. The geodesic congruence corresponding to 
such a universe includes geodesies with velocities arbitrarily close to the speed of light relative 
to C. Fast moving observers in this congruence would initially detect an extremely high rate 
of bubble hits. The short period of inflation inside the bubble slows down the congruence 
to non-relativistic speeds, and the rate approaches a constant on the hyperboloid (at large 
distances from the origin). Nevertheless, the angular asymmetry in the arrival directions 
of the bubbles remains. The reason is that, to an observer far from the origin, the initial 
surface looks very anisotropic, sloping down in time in the radial direction further away from 
the origin. In the unlikely event that we detect a signature of a bubble collision in some 
direction in the sky, we will be able to say that we are probably moving in that general 
direction relative to the preferred congruence C. 

Now that we know that the initial conditions at the beginning of inflation have a lasting 
effect on the distribution of bubble hits, one cannot help wondering what other effects they 
may have. We hope to return to this issue in the future. 
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Appendix 1 

In this appendix, we repeat the calculation of the distribution of collisions which we 
presented in Section Ha, but now instead of looking around the origin we consider the 
vicinity of a point £ D b s 3> 1. 

It is useful to go to the boosted frame where the point of observation is at the origin, as we 
did in Section II. Furthermore, we may take the limit 7 = cosh£ b s — > 00, (3 = tanh£ b s — > 1 
in the expression for the initial surface (14T1) . 

, _ cose' + F' 

ai ~ 1 + 2F'cos#' + F' 2 ' 

We may take Eq. ()42l) as the starting point. Performing the integral over a' we have 

3 / 1 77/2 \ 3 



dV A 



3F> 



1 + F' 
1 -F> 



l-F'- 



1 + 2F'cos#' + F' 2 



dF'dD! 



(78) 



Integrating over solid angle, we have 



^= fa(3+ r. +3f !.W. 



3(1-F') 3 (1 + F') 



(79) 



Finally, integrating over F' we obtain 



\4(e',^,^^-oo,Us>i; 



All 



IF' 1 - F' 

— m 

1 + 



3(1 — F' 



(80) 



After some algebra, this greatly simplifies to 



V^',t,Us » 1) = y [2e«'sinhe'tanh 2 (r/2) +e' 

This can be compared with the distribution around the privileged point £ Q b s = given in 
(121]) . The expressions are rather similar. In particular the asymptotic expressions (r ^> 1) 



(81) 
(82) 



are the same. 
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Suppose we match the inflating phase to a standard cosmological phase with scale factor 
given by a(r). As the universe decelerates, the horizon becomes larger and more bubbles 
come into sight. We can then ask what is the rate at which bubbles become visible per unit 
time to the observer at £' = 0. This is proportional to the nucleation rate A, times the rate 
at which the relevant 4- volume enters the backward light cone from the point of observation: 



dN{r) 
dr 



A 



dV* 



dm 

dr 



(83) 



Here r e is the time at the end of inflation, and the last factor gives the rate at which the 
observable distance changes as we look back to the surface where inflation ends: 

dr 



a(r) 



For instance, if inside the bubble the vacuum energy is zero, and we go immediately into 
the curvature dominated regime a(r) = r, then the rate at which we would see new bubbles 
entering our horizon per unit proper time is given by 



dN(r) 



AttX 



(84) 



dr 3 

Since this grows without bound, we would be guaranteed to see some bubbles sooner or later 
if we lived in this Minkowski vacuum. 



Appendix 2 

Here, we explore the possibility of setting up initial conditions which do not break the 
residual 0(3, 1) symmetry of the bubble. This is a formal exercise whose practical utility is 
unclear. Bubbles nucleate at random points, and it is not possible to set up initial conditions 
which preserve 0(3, 1) symmetry for all of the bubbles in the ensemble. Nevertheless, the 
attempt may be illustrative. 

The simplest thing we can do in order to have 0(3, 1) invariant initial conditions, is to 
require that there are no bubbles inside the backward light-cone from the antipodal point 
A of the nucleation event N (see Fig. [2]). We are also assuming that N is in false vacuum, 
so the surface of initial conditions is the disjoint union of the backward light cones from A 
and N. The relevant four-volume of our interest is represented by the shaded area in Fig. [21 
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FIG. 2: Conformal diagram of a bubble in a de Sitter space. The relevant 4-volume to the past of 
a sphere of radius £ around the origin is shaded in grey. 



It is convenient to use the conformal closed chart, in which the metric reads 

ds 2 = — l — (-drf + dx 2 + sin 2 x^ 2 )- 

COS Tj 



The 4-volume is given by 



Pu.ri Ltt / S -^d X dr ] = 87T 

COS 4 7] 



du 



tt/2 



dv 



sin [v — u\ 



(85) 



(86) 



Iq COS 4 (!> + «) 

where we have introduced the change of coordinates n = v + u and x = v ~ u - Past light 
cones are labeled by u = const values. Performing the integrals, we have 



S7T 



T 



2 cos 2u 
sin 2 2u 



— In tan u 



(87) 



Finally, we must relate the values of u to the values of £ on the r = const, hyperboloid. 
Alternatively, we may relate them to the by now familiar variable F. 
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On the bubble cone we have x = 71 ~~ Vi an d hence the physical radius of two-spheres is 
given by 

sin x 

R = = tan r] = — cot u. (88) 

cos 77 

On the other hand, in terms of a and F coordinates we have 

R = {l + aF) = ^j, (89) 

where we have used that on the bubble cone a — 2/(1 — F) [see Eq. (1151) ]. Equating both 
expressions for R we have 

F -I 

tanu= FTT- (90) 



Substituting in (IHTj) we obtain 



v - 8n 
T 



4F(F 2 + 1) _ 2ln ^~F 



(91) 



_ (F 2 — l) 2 1 + F_ 

Note that the leading dependence for large £ is the same as in V4, given in (1201) . However, 
both expressions differ in the subleading terms, and so the limit of large £ b s studied in 
the previous Section does not exactly agree with the result which we have just obtained by 
setting up Lorentz invariant initial conditions. 

Formally, we have obtained a finite result which does not break the homogeneity and 
isotropy of the reference bubble, by assuming that the point of observation is still in false 
vacuum. However, it is clear that the 4-volume V4 diverges when we set r = in the lower 
limit of integration . What this means is that the bubble has probability 1 of being hit 
by other bubbles immediately after its formation. In other words, there is no possibility of 
eternal inflation with these initial conditions: we have included too much of the contracting 
part of de Sitter space. 

Finally, we could have adopted yet another approach which formally does not require 
any "cut-off" initial surface. Indeed, we could ask the following question. Given that a 
particular point has not been hit by a bubble at some late time in the true vacuum (e.g. 
by the time of last scattering), what is the expected distance to the nearest bubble? If the 
true vacuum is of sufficiently low energy, the relevant 4-volume is finite, without the need 
of any cut-off. This is illustrated in Fig. 3. For clarity, we draw the case where the "true" 
vacuum inside the bubble is Minkowski. In practice, it is enough that the true vacuum be 
of sufficiently low energy density, so that the observer's cosmological horizon is much bigger 
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than the Hubble size during inflation. In this case, the shaded area corresponding to the 
4- volume available for the nucleation of bubbles which will hit a distance A£ away from the 
point of observation, is finite. Formally, we avoid the need to specify the initial surface. 

Again, this formal set-up is unphysical. Some initial conditions must be specified, and 
these will look different from the point of view of different observers. This will lead to the 
anisotropies which we have discussed in the present paper. 
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FIG. 3: Same as in Fig. [21 but now the "true" vacuum inside the bubble is of lower energy. For 
simplicity, we take it to be Minkowski, although this is not essential. In this case, the shaded area 
corresponding to the 4- volume available for the nucleation of bubbles which will hit a distance A£ 
away from the point of observation, is finite. All we need is that the backward light cone from 
the point of observation reaches the point \ = at some time in the past. This will happen if 
the cosmological horizon at the time of observation is many times larger than the Hubble size 
during inflation. Formally, we do not need a cut-off initial surface, and seemingly we obtain a 
finite answer which respects 0(3, 1) invariance. However, we know that physically some initial 
conditions are needed. These look different to different observers in the FRW congruence, which 
leads to anisotropics in the distribution of bubbles. 
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